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Systems with strong spin-orbit coupling, which competes with other interactions and energy 
scales, offer a fertile playground to explore new correlated phases of matter. Weyl semimetals are an 
example where the phenomenon leads to a low energy effective theory in terms of massless linearly 
dispersing fermions in three dimensions. In the absence of interactions chirality is a conserved 
quantum number, protecting the semi-metallic physics against perturbations that are translationally 
invariant. In this letter we show that the interplay between interaction and topology yields a novel 
chiral excitonic insulator. The state is characterized by a complex vectorial order parameter leading 
to a gapping out of the Weyl nodes. A striking feature is that it is ferromagnetic, with the phase of 
the order parameter determining the direction of the induced magnetic moment. 
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The emergence of materials whose properties are 
strongly influenced by spin orbit coupling is an exciting 
new phenomenon in condensed matter physics. Topolog- 
ical insulators [H 0] and Weyl semi-metals are two 
canonical examples, where the degrees of freedom, de- 
scribing the low energy physics of a non-relativistic many 
body system, are linearly dispersing massless fermions. 
An interesting possibility in these materials is the emer- 
gence of new phases of matter Q. In this letter we fo- 
cus on the Charge Density Wave (CDW) and Excitonic 
Insulating (EI) states in Weyl semi-metals. 

Weyl fermions are conjectured to be the low en- 
ergy excitations of Pyrochlore Irradiates (PIs)[l] and 
Topological-Normal insulator (TNI) heterostructuresjl, 
0] . The semi-metallic nature is derived from the touching 
of two non degenerate bands at an even number of Weyl 
points. The existence of such points leads to a number of 
anomalies in transport characteristics [8l-[lo| . While the 
physics of Weyl fermions has been extensively studied in 
the context of liquid "^He Hl, ll^, where they arise in the 



A phase, the recent developments have renewed searches 
for other systems that realize the phenomena [is}. 

An important property of Weyl points is that they are 
robust to most interactions. This is a consequence of 
chirality being a conserved quantum number in three di- 
mensions for massless fermions. Only interactions that 
couple fermions of opposite chirality can open a gap. For 
systems that preserve inversions but break time reversal, 
every Weyl node, if it exists, is accompanied by another 
of opposite chirality in the Brillouin zone. In the case of 
Pis 24 nodes at symmetry related points are conjectured 
to exist [1] while the TNI heterostructures have two such 
nodesjl]. In both cases there is perfect nesting. Conse- 
quently we seek to address the question of the possible 
particle-hole instabilities promoted to by repulsive inter- 
actions. 

To highlight the physics of the novel new state, we 
simplify to the case of two Weyl nodes and local den- 
sity density interactions. A more general analysis includ- 
ing the effects of long range Coulomb will be reported 
elsewhere. There are two types of particle-hole excita- 



tions that can arise in this case i) intra-nodal (occurring 
at zero momentum) and ii) inter-nodal (occurring at a 
finite fixed momentum associated with the nesting vec- 
tor). These are excitonic phases[14], the former being 
the EI while the latter is the CDW. Unlike conventional 
condensates in these sectors, the electron hole pairing of 
Weyl fermions leads to chiral phases. A minimum in- 
teraction strength is required to nucleate them which is 
the consequence of the vanishing density of states at the 
node. For local interactions, the chiral EI has the low- 
est threshold, opens a gap at the nodes and is the most 
stable state. The sign of the gap is opposite at the two 
nodes, preserving inversion symmetry. 

The order parameter breaks SU{2) x U{1) symmetry 
corresponding to rotational and phase invariance. This 
is in addition to the time reversal symmetry broken by 
the parent state. To characterize the broken symmetry 
state we define an orthonormal basis for three dimen- 
sional space {l,rh,n\. In particular, for the cartesian 
coordinate system with z as the quantization axis, the 

ordered state is one with 



= A 



(^)exp(^x), 



is the fermion anni 



where c-'^ 



hilation operator at the two Weyl nodes, labeled L 
and R, at momentum k in the band labeled ±, e| = 



{—ky/^^kl + k^^kx/yk"^ + k^,0}, A is the magnitude 
and X is the phase of the order parameter. The loss 
of rotational invariance is reflected in the the vector 
e| appearing in the order parameter. In general it is 

a vector lying in the plane spanned by I and m (i.e. 



^1 



hxk/ 



X k 



with n = Ix rh) and the corresponding 



order parameter is A 



V2 



exp(zx). 



The order parameter does not break inversion symme- 
try. Therefore the spin orientation for electrons at mo- 
mentum k and —k are identical. The effect of the sym- 
metry breaking is to cant spins in momentum space along 
a direction determined by the phase of the order param- 
eter. Thus vortex lines, if stabilized, in this system have 
an associated spin texture with finite divergence, and the 
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FIG. 1: The interaction shown in Eq.Q is a function of three 
vectors (g, and e|) that form a right handed coordinate sys- 
tem. Each vector couples to an operator of distinct symmetry 
in the particle hole channel. 

singularity localized in the vortex core. 
Model : Consider a system with two Weyl nodes at Ki — 
Kqx (labeled R) and —Ki — —Kqx (labeled L) with 
chiralities +1 and —1 respectively. The Hamiltonian is 

Ho± = ±hv ^Ijap ■ (k T Koj ^kfi (1) 

k 

where v is the fermi velocity and a — {ax,cFy,az} is a 
vector of Pauli matrices. The dispersion at each node is 

Cq- = zLhv \ q\ centered around ii^o, with q — (^k ^ Kq^ . 

The conduction (valence) band at the R node has its spin 
parallel (anti-parallel) to while the opposite is true 
at the L node. The general particle particle interaction 
takes the form 

y ^Y.j drd7'V{r-r')il^l,{7')i^,,{T')i:l{r)iljAT) 
k.k' ,q 



For the moment we do not make any assumptions on 
the nature of the interactions. Since the Weyl physics 
is the low energy description of a more general theory, 
we enforce an upper cutoff in the momentum integrals 
(up to an energy A) around the Weyl point. The precise 
source of the interaction and the renormalization effects 
will be considered elsewhere. Here we focus on the 
instabilities within a mean field analysis. 

Particle Hole instabilities : We rewrite the interaction 
in the basis of the non-interacting bands. Define ip^'^ = 

V^±c^'± as the fermonic field of the four bands, where -q 
is the spinor and c is the fermion annihilation operator. 
Of the 16 possible terms from Eq.([5]) only 6 terms satisfy 
momentum conservation for scattering restricted to the 
states within the cutoff around the node. 

For every momentum q = qq, where q = {qxi %, qz} is 
the unit vector along q, we define two orthogonal vectors 

4 = (>q = {^xqzl ^ql + ql, qyqz/^ql + ql,-^ql + ql} 

and el = 0,- = {-%! ^ql + qlAxI ^ql + g^'O}' 
that g, ei and e?. form a right handed coordinate system 
(see Figll]). The unit sphere is spanned by the vector q by 
two rotations, one about any axis perpendicular to e|.and 
the another about e?.. For example if we choose the first 
to be the z-axis, than the vector e?., which is the is in 
the spherical polar system, spans a unit circle (perimeter 
of the shaded region in FigH]) and the vector ei, which is 

the corresponding 0, spans the southern hemisphere. The 
following construction holds for an arbitrary quantization 
axis n, with the corresponding polar and azimuthal angle 
for q defined in the coordinate frame {Z, m, n\. In the rest 
of the letter we use the {i, y, z} coordinate system. 

Specializing to potentials that are even functions of fc, 
i.e. V^(fc) = V^(— fc), the interaction in terms of eg = 
el -f ie\, is 
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(3) 



The first and second term promote CDW in- 
stabilities with intra-nodal order parameter, i.e. 
i^tArc'tc'^ ^) 7^ with Ar an odd function of 

\ nfc ~^nfC I 

k, while the last term leads to inter-nodal order with 

Inter-nodal Charge density wave : We begin by study- 
ing the inter nodal instability that establishes ordering 



at 2Ko (third term in Eq.([3])). For momentum indepen- 
dent interaction potentials, V{k) = g/fl, where il is the 
volume of the system, the coupling takes the form 

^// = -iEE(^^:^"„.)-(^'^"I.^^.) (4) 

k,k' "=± 

Eq.Q is identical to that of the interaction in '^He 
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in the particle particle channel that leads to chi- 
ral superfluidity [ill [l^ . In Weyl semi- metals the 
corresponding state in the particle-hole channel is a 
charge density wave. Within mean field there are 
two possible instabilities: i) Chiral CDW: Ac = 

fi(EffHy-nk) = ^^(^) a'ld ii) Polar CDW: 



— ApZ. Note that the direc- 



P n \L^k' ~nk'~-nk' 

tions chosen for the ground state is for convenience and 
all other that are obtained by rotation in three dimen- 
sions are equivalent. The former is a chiral state while 
the latter is a non-chiral p wave CDW. The general form 
of mean held Hamiltonian is 
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FIG. 2: Gap magnitude and the change in free energy of 
the CDW. Left; Order parameter magnitude as a function of 
inverse of the interaction strength for chiral (blue o) and polar 
(purple □) states. Right: Comparison of condensate energy 
Ec with /i = for aforementioned states as a function of the 
interaction strength. 
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with A representing Ac or Ap. Minimizing the free en- 
ergy with respect to the order parameter, the integral 
equation that determines the condition on the coupling 
constant leading to the instability is 



lA 



c/p 



(5) 



Here Er 

k 



{hv\k\)^ + |A • fcp, /3 = l/ksT, Ac/p rep- 
resents the unit vector direction of chiral or polar state, 
and Vl = (■^) ; where n is the number of electrons 
in the system, set to 1 in the following. The summation 
over momentum is cutoff by the scale A/hv. In the limit 
of A ^ A and zero temperature, to leading order Eq.([5]) 
simplifies to 
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Both equations yield the same critical value for the cou- 
pling constant. For g > 3{hv)^ /2TrA'^ there exists charge 
density wave instability within mean field. Crucially the 
CDW instability does not gap out the Weyl nodes. 

In Fig[2] the zero temperature free energy difference 
between condensate and normal states, denoted as Ec or 
condensation energy in the figure, and magnitude of the 
order parameter of the two states are plotted as a func- 
tion of the interaction strength. As noted analytically, 
the two have the same instability threshold. The magni- 
tude of the order parameter is larger for the chiral state 
as compared to the polar state for the same interaction 
strength, leading to a greater (more negative) condensa- 
tion energy for the former. Thus in this sector the lowest 
energy state is the chiral CDW. 



The same conclusion is reached by studying the 
Ginzburg Landau (GL) theory, for a single order param- 
eter, near Tc. For T ~ Tc the free energy difference is 
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For the chiral CDW phase the coefficients are 
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The coefficient a in Eq.([7|) is to leading order in (/?cA)~^ 
and '^{x) is the digamma function. The critical tempera- 
ture Tc, evaluated by taking A ^ at T = Tc in Eq.dSJ, 

is ksTc c 



.(A2 



For the polar CDW phase 




the coefficient a is the same and b is 3/2 times the cor- 
responding coefficients in Eq.([7]). Thus the critical tem- 
perature Tc is the same as the chiral CDW phase. The 
coefficient of the quartic term is smaller for the chiral 
phase and hence it has the larger order parameter and 
lower free energies. The most stable state in the inter 
nodal sector is the chiral charge density wave. 
Intra-nodal Excitonic Insulator : Having established the 
instabilities in the inter-nodal sector we turn to those 
promoted by the first two terms in Eq.(|3]). For V{q) = ^ 
we get 
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There are four possible particle hole instabilities for 

i) Chiral-z EI: 

^. )> = 



order parameters with ei component 



(A^,i + A«i)(^), ii) Polar-z EI: (A^.^ -f A^,,)z, 
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FIG. 3: Left: Order parameter magnitude as a function of 
inverse of the interaction strength for Polar-x EI Ap^i (blue 
o), Chiral-z EI A^zi (purple □), Chiral-x EI Acxi (brown 
O), and Polar-z EI Ap^i (green A). Right: Comparison of 
condensate energy Ec with fi = for aforementioned states 
as a function of the interaction strength. 
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FIG. 4: Left: Order parameter magnitude as a function of 
inverse of the interaction strength for Polar EI Ap2 (blue o), 
Chiral CDW Ac (purple □), and Chiral EI Ac2 (brown o). 
Chiral CDW is shown here for comparison. Right: Compar- 
ison of condensate energy Ec with /i = for aforementioned 
channels as a function of the interaction strength. Chiral EI 
is the most energetically favorable among all states. 



iii) Polar-x EI: (A^^^ + ^pxi)^^ and iv) Chiral-x EI: 
(A^^i + Af^i)(^). This is related to the fact that el 
spans only the southern hemisphere and makes the z dif- 
ferent from the other two axes. The zero temperature or- 
der parameter magnitude, A^ — 2Aq^^, where a denotes 
the four possible phases, is determined by the same min- 
imization condition as in Eq.(l5]). For g > 3{hv)^ /2ttA'^ 
there exists polar-z EI which has the largest gap for the 
same interaction strength among these four states. At 
zero temperature we get Ap^i > A^xi > A^^i > Ap^-i 
and EP'^^ < Ef"^ < E^'-^ < E^^'^ for the same inter- 
action strength g larger than 6(ftw)^/7rA^ as shown in 
FiglSl By comparing it with Figl2]we see the polar-z EI 
and the chiral CDW are equally energetically favorable 
among the six possible states aforementioned for a given 
interaction strength. The polar-z EI state has Ginzburg 
Landau parameter a and 6 as i and j of the chiral CDW 
state, giving rise to smaller order parameter {1/^/2 that 
of the chiral CDW) but the same free energy as chiral 
CDW at finite temperature. None of these phases gap 
out the Weyl node. 

For order parameter of particle hole instabilities along 
e| component the order parameter have different signs 
between the two Dirac nodes and there are two possi- 
ble EI states based on symmetry (no z-component in 
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)) = (A^2 + ^c2)i^) with the signs fixed by 
setting A^ = A^2 = ^c2/2 > under the assumption 
of inversion symmetry, and ii) Polar EI: (A^2 + ^a2){^) 
with the same sign convention. At zero temperature, the 
chiral EI gap equation Eq.([5]) yields the condition 
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For g > {hv)^/TrA'^ there exists a chiral EI instability 
with a uniform gap. From Fig|4] we see this state has 
the largest gap value compared with all other state for 
a given interaction strength. The chiral EI state is also 
more energetically favorable at zero temperature com- 
pared with the chiral CDW or polar-z EI. It is the only 
state that opens a gap at the Weyl nodes. 

The chiral EI states mix particle and hole states which 
have opposite spin orientations in the non-interacting 
limit. Thus the ground state no longer has the spins 
aligned with their momenta. To evaluate the nature of 
the spin configuration we compute the expectation value 
of spin at momenta k for the occupied band. The result 
is that there exists a net polarization at each Weyl node. 
To understand the origin of this effect we first rotate the 
mean field Hamiltonian back to the ^'k.a basis. It takes 
the form 
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where A' = A sin (x) I + A cos (x) to, A is | Ac2 ' ^^1, and 
X is the corresponding phase. Under inversion e| does 
not change sign but e| does. Since we also go from one 
Weyl node to the other under inversion, the Hamiltonian 
preserves the symmetry. Averaging over polar and az- 

imuthal angle the Hamiltonian is A • (t/2 at both nodes. 
This is the origin of the magnetization in the chiral EI 
state and serves as a diagnostic of the state. An external 
magnetic field couples linearly to the order parameter 
leading to a lowering of the critical coupling, provided 
that it is smaller than magnitude of the gap and no Lan- 
dau levels are formed. 

Discussion : The interplay of strong spin-orbit coupling 
and repulsive interaction leads to novel excitonic phases 
in Weyl semi-metals. For the simple case of two Weyl 
nodes there are both intra-nodal and inter-nodal insta- 
bilities that can be realized. For short range interaction 
the fully gapped chiral excitonic insulator has the lowest 
free energy. The p-wave nature leads to vectorial order 
parameters. Crucially one can define three orthonormal 
vectors each of which couples to three distinct order pa- 
rameters namely CDW, inversion preserving EI and in- 
version violating EI. In analogy with superfluid liquid He, 
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the finite temperature phase diagram will be determined 
by the interplay between the various sectors, which has 
been ignored our GL theory for single order parameters. 
In addition the phase of the EI order parameter deter- 
mines the direction of magnetization at each node. Thus 
vortices have a magnetic character resulting in spin tex- 
tures. The precise nature of these excitations and nature 
of the phase diagram, both for local interactions and long 
range Coulomb interactions, is currently under investiga- 
tion. In particular, in analogy with '^He A, it is of great 
interest to determine if any of the nodal excitonic phases 
are possible at finite temperatures. 

Our analysis focussed on systems with two Weyl nodes. 
The TNI heterostructures[^] offer a simple realization 
where one can look for the excitonic phases reported here. 
Pis are another system where, in a regime of intermediate 
correlations, a Weyl semi-metal with 24 nodes is conjec- 
tured to occur 0. This implies that each Weyl node can 
couple to 12 other nodes with opposite chirality. On the 
other hand the available phase space and the cut off A 
will scale down. Whether the parameters end up being 
favorable to obtain the state is an open question that 
can only be answered once a Weyl semi-metallic state is 
established and characterized. Nevertheless the on going 
efforts to investigate systems with strong spin orbit in- 
teractions promises to focus interest in exploring possible 
new correlated phases of matter. 
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